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We study possible paramagnetic phases of antiferromagnets on the Shastry-Sutherland lattice by 
a gauge-theoretic analysis of fluctuations in a theory with Sp(2iV) symmetry. In addition to the 
familiar dimer phase, we find a confining phase with plaquette order, and a topologically ordered 
phase with deconfined S = 1/2 spinons and helical spin correlations. The deconfined phase is 
contiguous to the dimer phase, and in a regime of couplings close to those found in the insulator 
SrCu2(B03)2- We suggest that a superconductor obtained by doping this insulator with mobile 
charge carriers will be an attractive candidate for observing the anomalous magnetic flux properties 
associated with topological order. 
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I. INTRODUCTION 

Much interest has recently focused on-the magnetic 
properties of the insulator SrCu2(B03) j!ra. The low 
energy spin excitations in this material reside on the 
S = 1/2 Cu ions which lie in two-dimensional layers 
decoupled from each other. The experiments show a 
clear indication of an energy gap towards spin excita- 
tions, making this one of the few known two-dimensional 
systems with a spin gap. Remarkably, the pattern of 
near-neighbor antiferromagnetic exchange couplings be- 
tween the Cu ions turns out to be identical to that in 
a model Hamiltonip|n studied many years ago by Shas- 
try and SutherlancH. These authors also showed that a 
simple decoupled dimer wavefunction was an exact eigen- 
state of this Hamiltonian, and that it was the ground 
state over a restricted parameter regime. 

The Shastry-Sutherland antiferromagnet is sketched in 
Fig [jj The Hamiltonian is 
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where S, are S — 1/2 operators on the sites, i, of a 
square lattice. The exchange J\ > acts along the near- 
est neighbor links (shown as full lines in Fig |l|), while 
J2 > acts on the diagonal links, shown as dashed in 
lines in Fig [l]. It was established!] that a simple prod- 
uct of singlet pairs on the diagonal links was the ground 
state of H for sufficiently large 3i/J\. However, an un- 
derstanding of the experiments requires a description of 
the excitation spectrum, and also of possible quantum 
phase transitions to other states at smaller J2/J1. These 
issues have bees, addressed in a numbeipefrjecent theo- 
retical worksoiii Many of these studiesQETH involve nu- 
merical analyses based upon large-order series expansions 
departing from various decoupled cluster states. Quan- 
tum Monte Carlo simulations have, in principle, a smaller 



bias due to the choice of an initial state, and can be ex- 
tended to much large system sizes; however, simulations 
of H suffer from a sign problem, and so such studies have 
not been possible. An analytic mean-fielcLapproach has 
also been undertaken by Albrecht and Milafl: results were 
obtained mainly for the magnetically ordered states, and 
the various distinct paramagnetic states were not distin- 
guished. 

Quite apart from determining the ground states of the 
specific Hamiltonian H, it is also of interest to determine 
the phases of models which are "near" the parameter 
space of H . This is in the hope that future experiments 
may succeed in deforming the insulator SrCu2 (1303)2 by 
substitutional doping (which can induce mobile carriers), 
or by the application of hydrostatic pressure. Doping the 
antiferromagnetic insulator La2Cu04 led to the discovery 
of high temperature superconductivity: related phenom- 
ena may be expected here, although, as we shall argue 
later, the presence of strong frustration in the parent in- 
sulator SrCu2(BOa)2 may lead to profound differences in 
the nature of a possible superconducting ground state. 

This paper will examine a generalization of H to 
Sp(2A) symmetry (SU(2) = Sp(2)) and describe the 
properties of the large N limit. Some of the phases 
obtained in such a large N limit may not actually ap- 
pear in the phase diagram of the SU(2) model H — nev- 
ertheless, as we have just argued, the phases may still 
be of relevance to physical systems whose microscopic 
Hamiltonians are near the parameter space of H. Such 
an approach has been fruitfully applied to a number of 
other_fmstrated quantum antiferromagnets in previous 
workETtLS. The method leads to an unbiased selection 
of possible ground states in the large N limit, both with 
and without broken spin rotation symmetry. Moreover, 
a gauge-theoretic description of the fluctuations about 
the mean-field solution allows a systematic and reliable 
assessment of the stability of the various ground states, 
along with a description of the dynamics of the excita- 
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tions. 

The Sp(27V) generalization of H is defined by introduc- 
ing canonical Bose creation operators b\ a on every site i, 
with a = 1 . . . 2N a Sp(27V) index. The allowed states in 
the Hilbert space satisfy the constraint 



bLb? = 2NS 



(1.2) 



on every site i (we follow the convention of summing ove r 
all repeated Sp(2iV) indices); the right hand side of ( |l.2j ) 
must be a positive integer, the values of S are constrained 
accordingly — for the physical case, N = 1, S must take 
half-integral values, as expected. The Hamiltonian is 
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where J a P — J a p = —Jp a is the generalization of the 
antisymmetric e tensor of SU(2) (i.e.) J contains N 
copies of e along its center block diagonal, and vanishes 
elsewhere) . 

The large N analysis of a large class of models, of which 
H is a member, was described with some generality in 
Section II of Ref. [H]. We will follow the same method 
here, and so will dispense with the details of the compu- 
tation. The resulting mean-field phase diagram is shown 
in Fig U as a function of J 2 / ' J\ and 1/5 (in the large N 
limit, S becomes a continuous real variable). The posi- 
tions of the various phase boundaries are not expected to 
be quantitatively accurate for the physical N = 1 case. 
However, the general topology of the phase diagram, the 
nature of the phases and their excitations, and the crit- 
ical properties of the quantum phase transitions can be 
reliably described using Fig ^ as a starting point. 

The properties of all the phases in Fig || will be dis- 
cussed in detail in Section |J. Here we highlight our main 
new results. 

One of the paramagnetic phases has short-range equal- 
time spin correlations peaked at the wavevector (tt, 7t). 
[We denote this phase (tt, tt) short-range ordered (SRO) 
in Fig [| here we are placing the sites on the vertices 
of a regular square lattice as in Fig [l], and measuring 
wavevectors in units of 1/ (nearest neighbor spacing). In 
the experimental SrCu2 (603)2 system, the positions of 
the sites is different, and there will be a corresponding 
transformation in the wavevector dependence of observ- 
ables.] At the mcpf-Eeld level, this phase is identical to 
that found earlierEBEEl on the square lattice with J 2 = 0. 
However, we will show here that a difference does emerge 
upon caaaideration of fluctuations. For J 2 — 0, it was 
shownllacJ that Berry phases associated with hedgehog- 
instantons led to columnar spin-Peierls order in the (tt, tt) 
SRO phase. Here we show that a closely related anal- 
ysis for the Shastry-Sutherland lattice leads instead to 



"plaquette" order in this phase. Just such .a phase was 
considered recently by Koga and KawakamiQ. 

Our other new results are also associated with a param- 
agnetic phase. This phase is denoted (tt, q) SRO in Fig || 
and is obtained by a destroying the long-range magnetic 
order in a helically ordered phase. Equal-time spin corre- 
lations show short-range incommensurate order, and the 
spin structure factor is peaked at the incommensurate 
wavevector (tt, q) (the value of q varies continuously as a 
function of J2/J1). As in previous incommensuraJtc-SRO 
phases found on frustrated square lattice modelsE3'L3, we 
argue that the excitations above the ground state are de- 
confined spinons which cajry spin 5=1/2 (for SU(2)). 
Also as in previous workcallil, the quantum phase transi- 
tion between this phase and the (tt, tt) SRO (plaquette) 
phase (Fig |^) is described by theory of a charge 2 Higgs 
scalar coupled to a compact U(l) gauge field; the decon- 
finement transition is associated with the condensation 
of the Higgs field, aj^i^J^CEilM^ properties are those 
of a Z 2 gauge theory^'ll^nlJEJa. We will also consider 
here the transition between the deconfined phase and the 
dimer phase: by a somewhat different analysis, we will 
show that this transition also reduces to a Z2 gauge the- 
ory description. 



II. MEAN FIELD PHASE DIAGRAM 



As discussed in Section II of Ref. 15, a key quantity 
determining the nature of the phases is a complex, di- 
rected, link field Qij = —Qji. Operationally, this field 
is introduced to decouple the quartic boson interactions 
in H by a Hubbard-Stratonovich transformation. After 
this decoupling, the effective action contains the terms 



S= dr 
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where r is imaginary time, Jy = J\ (Jij = J2) on the 
horizontal/ vertical (diagonal) links, and the ellipses rep- 
resent standard terms which impose the c anon ical boson 
commutation relations and the constraint (1.2). It is also 
clear from the structure of S that the average value of 
On satisfies 



(Qv) = ±{j aP bLbl 
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(2.2) 



For larger values of S, the dynamics of S requires con- 
densation of the bf bosons, and hence a non-zero value 
of 



(2.3) 



such phases break the spin rotation symmetry, and have 
magnetic long-range order. As described in Ref |l5|, we 
optimized the ground state energy with respect to varia- 
tions in (Qij) and xf for different values of J2/J1 and S. 
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The four-site unit cell of the Shastry-Sutherland lattice, 
depicted in Fig |[ has 10 different Qij fields. Care must 
be taken to identify gauge-equivalent configurations. We 
find that each saddle point may be described by purely 
real (Qij}- The resulting phase diagram is shown in Fig|jj. 
We describe the phases in turn in the following subsec- 
tions, considering first the mag netically ordered phases 
with xf 7^ in Subsection II A, and then the paramag- 



netic phases in Subsection II B 



A. Magnetically ordered phases 

1. Neel (it, it) LRO state 

This is the familiar long-range ordered (LRO) state in 
which (Si) is collinearly polarized in opposite directions 
on two checkerboard sublattices. It is known to be the 
ground state of H for J 2 = 0, S = 1/2 in the physi- 
cal N = 1 limit. A gauge may be chosen in which the 
expectation values of link variables, (Qij), are nonzero 
and equal on the horizontal and vertical links, while the 
expectation values on the diagonal links are zero. In 
the notation of Fig || then Qi — Pi (i — 1,2,3,4) and 
Ri = R 2 = 0. 



2. Helical (n,q) and (q,n) LRO states 

This phase is characterized by non-zero values of (Qij) 
on the horizontal, vertical and diagonal links. A gauge 
choice sets all the Qi equal to each other, and similarly 
for the Pi ; in the appendix we present an argument which 
shows the values of the Pi and Qi are also equal to each 
other. There are two gauge non-equivalent choices for 
the values of Ri t 2- one state has R\ = R2 and the other 
Ri = — i?2- The two states are interchanged under a 90° 
rotation, and correspond to spirals ordered in the hori- 
zontal or vertical directions. At large values of the spin, 
this phase appears at J2 > J\ , in accord with the classi- 
cal calculation of Shastry and Sutherland!! Equal-time 
spin correlations exhibit long-range incommensurate or- 
der, and the spin structure factor peaks at the incom- 
mensurate wavevectors (n, q) or (q, it) with the value of 
q varying continuously as a function of J 2 / Ji ■ This state 
also appears in the studies of Refs |£jjic|. 



B. Paramagnetic phases 

In this subsection we discuss the three phases for which 
xf = 0. As a consequence, spin rotation symmetry is 
preserved and only spin SRO arises; however there may 
be ordering in other singlet order parameters. 



1. Dimer state 

This is the exact SU(2) eigenstate of decoupled singlet 
pairs found by Shastry and Sutherland! In the large 
TV limit, this corresponds to a saddle point at which the 
(Qij) are non-zero only on the diagonal links: R\ — R2 ^ 
and Qi — Pi = 0. Note that the bf boson are spatially 
decoupled at such a saddle point: each bf can only hop 
across a single diagonal link. This simplifies analysis of 
fluctuations about the saddle point in or near the dimer 
state, as will be discussed in Section IV. At higher orders 
in 1/N, the bf can indeed hop through the entire lattice; 
we expect that the lowest lying excitation will be a S = 1 
spin tripleta (for SU(2)), consisting of a confined pair of 
the bf bosons. 



2. (iT,1t) SRO 

This state is obta ined by quantum-disordering the Neel 
state of Subsection II A 1 , and the expectation va lues o f 
Qij have the same structure as those in Subsectio n |ll A 1 . 
As has been discussed in some detail in Refs the 
quantum fluctuations in this phase are described by a 
compact U(l) gauge theory. Such a theory is always 
confining, and thus the bf bosons again bind to yield 
a S = 1 quasiparticle above a spin gap. There is also 
an interesting structure in the spin-singlet sector: this 
is considered in Subsection III where it is demonstrated 
that at finite N this phase has "plaquette" order. 



3. (ir,q) and (q,ir) SRO 

In this phase (Qij } are non-zero on the diagonal, hori- 
zontal and ver tical lin ks, like the helical (it, q) LRO phase 
of Subsection |ll A2[ Again there are two gauge non- 
equivalent configurations, corresponding to the choices 
Ri = R2 with Qi = Q < Pi = P [the (ir,q) phase], 
and Ri = -R 2 with Q 4 = Q > Pi = P [the (q,n) 
phase] . Thus all of the horizontal Qi fields acquire the 
same expectation value, but unlike in the helical LRO 
phase, this value differs slightly from that of the vertical 
Pi fields; the difference is only on the order of a part in 
ten thousand. The state is a spin-singlet and there is a 
gap to all spin excitations. Nevertheless, the symmetry of 
90° rotations between the vertical and horizontal direc- 
tions is broken — this would now be apparent in various 
spin-singlet observables like the bond exchange energies 
or the bond-charge densities. This phase may therefore 
be viewed as a spin-singlet "nematic" as only rotational 
symmetry is broken. The choice of a vertical or horizontal 
spatial polarization in the nematic order leads to a two- 
fold degeneracy in ±h|e_mjmMr^ate. The state also has 
"topological" ordeiSBBHO'Ea, and this would lead to 
an additional four-fold degeneracy in a torus geometry. 
Unlike the commensurate SRO phases, the spinons are 
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deconfined. We describe the deconfinement transition 



below in Section IV. The spinon dispersion has its min- 
ima at momentum (ir/2,q/2) or (q/2,7r/2). Although 
this phase is realized only for S < 1/2 in the large N 
limit, it seems possible that in the physical limit N = 1 it 
could extend up to S = 1/2 for a narrow range of Jij J\- 
Similar behavior was found in a study of the Sp(2iV) 
Heisenber|g-.antiferromagnet on the anisotropic triangu- 
lar latticed. It would interesting to search for this phase 
using numerical methods. 

We conclude this section by briefly comparing our re- 
sults to other published calculations. For S = 1/2 we 
find that the transition between Neel and Helical LRO 
phases is continuous, occurring at J2/J1 « 1^2, close to 
the value of 1.1 found by Albrecht and MilaE, who also 
report a continuous transition. Also in agreement with 
Albrecht and Mila, we find that the transition between 
the Helical LRO and Dimer SRO phase is first order, but 
occurs sX J2/J1 ~ 2.7 instead of 1.65. Carpentier and 
BalcntgI9 also found a helical LRO phase, but presented 
arguments that an intermediate phase may exist between 
the helical LRO and dimer phases: our (tt, g),-SRO state 
is precisely such a phase. Koga and KawakamiS employed 
a series expansion to find, for S = 1/2, a plaquette phase 
which intervenes between the Neel and dimer phases. As 
shown below, the (tt, tt) SRO phase acquires plaquette 
order at finite N, but as can be seen in Fig 0, at large 
N this phase only occurs for S < 1/5. If finite N fluctu- 
ations push the phase boundary up to S = 1/2 then the 
following sequence of phases would occur as J2/ ( J\ + J2) 
increases from to 1: Neel, Plaquette (7r,7r) SRO, (tt, q) 
SRO, and finally Dimer SRO. 



III. PLAQUETTE ORDER IN THE 
COMMENSURATE PARAMAGNET 

This section will discuss the fate of the spin singlet sec- 
tor upon including fluctuations about the mean-field in 
the (tt, tt) SRO state. The results below are a straight- 
forward generalization of those obtained in Refs |l9|,^0] tor 
the square lattice antiferromagnet. We will only consider 
the case where 2SN is an odd integer (for the physical 
SU(2) case, this means that S is half an odd integer); 
the generalization to other values of S follows as in ear- 
lier work. 

In the present large N approach, regular perturba- 
tive corrections order by order in 1/iV do not qualita- 
tively modify the nature of the mea n-fie ld ground state. 
However, singular effects do appeart3E3 upon consider- 
ing the consequences of 'hedgehog' like instanton tunnel- 
ing events and their Berry phases. Such a calculation is 
technically involved, and a somewhat more transparent 
discussion of essentially the same physics emerges from 
studying the "quantum dimer" modelEj (see Appendix 
A of Rcf for a discussion of the equivalence between 
the instanton physics of the large N expansion and dual 



representations of the quantum dimer model). Here we 
shall follow the treatment of Ref |29|. 

The quantum dimer model represents the Hilbert space 
of low-lying singlet excitations by assuming that it can 
be mapped onto states represented by a near-neighbor 
singlet bond ('dimers') covering of the lattice. In the 
present (71", tt) SRO phase, we need only take dimers con- 
necting nearest neighbor sites on horizontal and vertical 
links. The dimers along the diagonal links are assumed 
to occur only rarely in this phase: they can therefore 
be integrated out, and serve mainly to modify the effec- 
tive Hamiltonian in the space of horizontal and vertical 
dimers. Indeed, the most important consequence of this 
procedure is apparent from a glance at Fig [l]: the diago- 
nal dimers divide the plaquettes of the square lattice into 
two classes, those with and without diagonal links across 
them, and we expect dimer resonance terms around these 
plaquettes to have distinct matrix elements (see Fig ||) . 
This distinction will be the only difference from earlier 
analysesliSEj, and we will show that it is sufficient to lead 
to plaquette order in the (tt, tt) SRO phase. 

Our results emerge from an analysis of , 
representation of the quantum dimer modeBEIffil 
There is a rigorous, one-to-one mapping between the set 
of coverings of the square lattice with nearest-neighbor 
horizontal and vertical dimers, and the configurations 
of an interface of heights, h a , defined on the sites, a, 
of the dual square lattice (we identify two interfaces is 
equivalent if they are related by a uniform translation 
ha — * ha + P, where p is any integer) . The values of h a 
are restricted to 



h a = n a + ( a 



(3.1) 



where n a is a integer which fluctuates from site to site, 
and Ca is a fixed fractional offset which takes the values 
0,1/4,1/2,3/4 on four dual sublattices, X,Y,Z,W, as 
shown in Fig S. We further restrict the h a to satisfy 
\h a — h \ < 1 for any pair of nearest-neighbor sites a, b. 
We can now specify the connection between the height 
model and the dimer coverings. Examine the value of 
\h a — h \ for every nearest neighbor pair, and if \h a — h \ > 
1/2, place a dimer on link shared by the plaquettes of 
the direct lattice around a and b. It is not difficult to 
see that a consequence of our choice of the £ a offsets is 
that dimers so obtained will form a close-packed covering 
of the lattice. Examples of the relationship between the 
height values and dimer coverings are shown in Fig |i[ 

We can now use general symmetry considerations to 
write down an effective action for the height degrees of 
freedom. As is standard in theories of interface mod- 



els, we promote discrete heights h a , in (3.1), to continu- 
ous real variables Xa by the Poisson summation formula, 
and "soften" the constraints to periodic cosine potentials 
whic h have minima at the values Xa — h a which obey 
( |3.l| ). In this manner we obtain the action 
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S x = dr 



(ab) 



K 

\ -Y( d TXa) 2 ~ Va COs(27T(Xa ~ Co)) 



. (3-2) 



where the sum over (ab) extends over nearest neighbor 
sites, and K is the stiffness towards spatial fluctuations of 
the interface height. The corresponding stiffness towards 
time-dependent fluctuations is K T , and, for simplicity, 
we have taken its value a independent. The symmetry 
of the lattice requires that the strength of the periodic 
potential take two possible values, y a = yi or y a — 7/2 de- 
pending upon whether the plaquette a has a diagonal J2 
link across it or not. This is the sole distinction from the 
analysis of the square lattice antiferromagnet in Ref |2C| , 
which had y\ — 2/2 ■ 

The fundamental property of interface models in 2+1 
dimensions, like S x , is that they are always in a smooth 
phase. This means that the symmetry of height trans- 
lations is always broken, and (x a ) = (h a ) has some def- 
inite value across the entire system. As was argued in 
Rcfs [l9],|o| any such definite value necessarily breaks the 
lattice symmetry of the underlying antiferromagnet, and 
will lead here to plaquette order. 

With the assumption of a smooth interface, the opti- 
mal interface configurations can be determine by a simple 
minimization of S x by a set of time-independent values 
of Xa- We allow for distinct expectation values, xw, Xx, 
Xy, and xz on the four dual sublattices. Then the prob- 
lem reduces to the minimization of the following energy 
as a function of these four real variables: 



£ v = K 



{Xx - Xw) 2 + (xw ~ Xy) 2 

+ (xy - Xz) 2 + (xz - Xx) 2 

- y t [cos(2ttxw0 - cos(2ttxf)] 

- m [sin(27rxx) - sin(27rxz)] 



(3.3) 



This minimization is a straightforward, but somewhat 
tedious, computation. The present analysis is valid only 
for small yi, j/2, and so we analytically determine the 
minima in power series in 2/1.2- We define 



XW = Xl + X2 + X3 
XX = Xi - X2 + X3 
XY = Xl + X2 ~ X3 
Xz = Xi - X2 - X3- 



(3.4) 



We find that at the saddle points of E x , 

^(yl 



X2 
X3 

Xa 



vi) 



sin(47rxO + 0(yf )2 ) 
sin(27r Xl ) + 0(yK) 



16K 2 
2K 

||cos(2 7 r Xl ) + 0(^ 2 ). 



(3.5) 



The average interface height, xi, is determined by the 
minimization of 



E x = E + Acos(4ttxi) 



Bcos(87rxi) 



(3.6) 



where Eq is an uninteresting constant independent of Xi , 



A = 



B = 



2K 

* 6 (7yt ^ 



6K 3 
7y 2 4 ) 
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(3.7) 



and all omitted terms are of order y\ 2 or higher (in ob- 
taining the results in (3.7) we had to include terms in 
( |3.5[ ) which are one order higher than those shown) . Note 
that the square lattice antiferromagnet, with y\ = y 2 , has 
A = 0. 

We now have to minimize (3.6) to determine xi- Then 
from (3.5) we know X 2,3,4i and hence the configuration 
of the interface heights. Then, from the connection be- 
tween I h a — hb I and the corresponding dimer occupation 
numbers, we can determine the pattern of the distribu- 
tion probabilities of the spin singlet bonds in the original 
antiferroma gne t . It is a simple exercise to determine the 
minima of (3.6) for different values of A and B; the re- 
sulting phase diagram is shown in Fig ^, and we now list 
the various minima and the associated ground states of 
the antiferromagnet. 

(i) A > 0, B < A /A: There are degenerate minima 
at xi — 1/4,3/4. The system spontaneously breaks a 
translational symmetry by choosing one of these minima. 
With the mappings above, it is easy to see that these are 
the plaquette states, one of which is depicted in Fig |6[ 

(ii) A < 0, B < — A/ A: Now the two equivalent min- 
ima are xi = 0, 1/2. These also correspond to plaquette 
states as above, but the chosen plaquettes are now around 
half of those containing diagonal links (see Fig ||) . 

(Hi) The remaining values of A and B have four degener- 
ate minima at xi = 1/4 ± z9, 3/4 ± d, where < 1? < 1/4 
varies continuously as a function of A/B. These states 
have spin-Peierls order of the type shown in Fig 0: the 
links are divided into four columnar sets, with each set 
having a different value of (S^ • Sj) on its links. This state 
interpolates between the plaquette state in (i) as d — > 
and that in (ii) as — > 1/4. 



The present analysis is for small j/i, and so, from (3.7) 
we should assume that B <C \ A\. Furthermore, the pres- 
ence of the frustrating J2 interaction on half the plaque- 
ttes means that the hedgehog tunneling events are more 
likely to be centered on these plaqu ettes . Using the map- 
ping of such ev ents to the model (3.2), we expect that 
yi > Z/2- From (3.7) we therefore conclude that the most 
likely possibility for the ground state is that in (i) above. 
The same state has also been considered in Ref |[ 

We conclude this section with a few comments on the 
(tt, 7r) SRO phase of the antiferromagnet with full square 
lattice symmetry, in which there is a diagonal J2 ex- 
change between every pair of next-nearest-neighbor sites. 
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Recent numerical work on such an antiferromagnetBEl 
has found evidence for spin-Peierls ordering with the 
same spatial structure as in (Hi) above for the Shastry- 
Suthcrland antiferromagnet. However, we noted earlier 
that the square lattice symmetry implies that A = 0: for 
this value, i? = 1/8, and the spin-Peierls state of (Hi) has 
a larger symmetry (two of the four sets of columnar links 
are equal to each other), and becomes equivalent to the 
ordering discussed in Refs ^9 20. To obtain $ ^ 1/8, and 
so a ground state with the symmetry of that in Fig ||, 
we need to add to E x a higher order term C cos(167rx): 
then there can be an eight-iold degenerate ground state, 
with d and 1/4 — d equivalent to each other. This is the 
state that appears to have been found in Refs |32] , |l8| . 

Note also that for the square lattice case, the B < 0, 
A = solution iias the four plaquette states degenerate 
with each otherEl 



IV. DECONFINEMENT TRANSITION OF THE 
DIMER PHASE 

The deconfined, "spin-liquid", (ir,q) SRO phase in 
Fig H is flanked on both sides by confining paramagnetic 
phases, the plaquette and the dimer phases. 

As we indicated Section || the deconfinement- 
confinement quantum phase transition from the (ir, q) 
SRO phase to the plaquette phase can be described in 
a theory essentially identical to that considered ram- 
ously for frustrated square lattice antiferromagnetsllBEJ. 
At the mean-field level, the transition is signaled by the 
onset of non-zero expectation values of Qij on the diago- 
nal links: we will denote these diagonal Qij as Qfj. Upon 
considering fluctuations, we find that the Qfj constitute 
a charge 2 Higgs field in a compact U(l) gauge theory, 
and the deconfi nmpn tp CQnJ lBement transition is that in a 
Z 2 gauge theoryyyiil&y. 

This section will consider the second deconfincmcnt- 
confinement transition in Fig |3|, between the dimer and 
(ir, q) SRO phases in more detail. We will see that this is 
also described by a Z 2 gauge theory, and the emergence of 
the Z 2 gauge symmetry can be described in a somewhat 
more transparent manner. 



As noted in Section II B 1 , the dimer phase is char- 
acterized by non-zero expectation values of the diagonal 
Qfj links. These links are all decoupled from each other, 
and this leads to a simple, local structure in the effec- 
tive action for the fluctuations. The transition to the 
deconfined phase is now signaled by the onset of non- 
zero expectation values of the Qij on the horizontal and 
vertical links, and we will denote these by Qfj and Qfj 
respectively. Near the phase boundary, we need only con- 
sider the structure of the effective action as a functional 
of the QjJ v after all other degrees of freedom have been 
integrated out. 

The simplest terms in the effective action arise from the 
on-site propagation of the bf on the site i in imaginary 



time. Integrating out the bf in powers of the Q^ 
lowest order terms have the form 



the 
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c 2 {Q12Q2Z Q^QS+H.c.} 

□ 



(4.1) 



where c\, c 2 are constants, the first sum is over nearest 
neighbor links, and the second sum is over plaquettes, 
with the sites labeled as in Fig [?]. A crucial property 
of Si is that all terms are invariant under a local U(l) 
gauge transformation 



' 'J 



t>i+4>j) 



(4.2) 



where the phase (f>i can take arbitrary distinct values on 
the sites i. 

We have so far not made use of the fact that the 
nonzero value of (Qfj) allows the bf bosons to hop across 
a single diagonal link. Such hopping processes will induce 
a large number of additional terms between the Qff. We 
will now write down the structure of all such terms which 
appear at fourth order in the Qfj 1 '. It is convenient to 
group these terms into sets associated with links emanat- 
ing from a given plaquette which does not have a diago- 
nal dimer across it: one such plaquette is that with the 
sites 1,2,3,4 in Fig ^ and we now write down all four- 
link terms in which every link has at least one site on the 
central plaquette. It is not difficult to see that all other 
four-link terms can be obtained by a simple translation 
of these terms to other plaquettes. The terms are 



<S 2 



ch 



+Q37Q84Q41Q12 — Q84Q51Q12Q32 + C.C. J- 

+C4IQ51Q26Q12Q41 — Q2QQ37Q32Q12 

+ Q37Q84Q43Q32 — Q 84Q 5lQ 4lQ 43 c - c - } 

— C5 ^Q^iQ 12 QstQ 43 + Q26Q32Q84Q41 + C.C.| 
+c 6 {Q h 51 Q V 2 6 Q h 37 Q v 8i + c.c.} . (4.3) 



Clearly, (4.3) is not invariant under (4.2). However, a 
residual Z 2 gauge symmetry does survive. We see that 
(4.1, 4.3), and all other allowed terms, are invariant under 



Q 



h , 1.' 



■\h,v 

hj ViVj 
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where rji — ±1 performs the gauge transformation. How- 
ever, it is not possible to choose the rji independently on 
every site: it is easy to see that we need the additional 
constraint 



G 



rji = rjj whenever i and j 

are separated by a diagonal link. 



(4.5) 



So the Z 2 gauge degree of freedom is halved from that 
present on the original square lattice. 

To place the Z 2 gauge theory in a more conventional 
form, it is useful to introduce a slightly different param- 
eterization of the degrees of freedom. First, we neglect 
all amplitude and phase fluctuations and replace all the 
Qij by discrete Ising variables taken only the values ±1. 
Then we choose to represent all the Qy as Ising gauge 
fields, a, while all the are written as products of a 
and a second Ising spin field, //; thus: 

Q" ~ tr/x. (4.6) 

This is shown a more explicitly in Fig (7]. Notice that each 
pair of horizontal and vertical links that form a triangle 
with a single diagonal link share the same Ising gauge 
field a. This choice is a consequence of the constraint 
(4.5) — as a result, all the fj, fields are invariant under 
the gauge transformation generated by the r]i, while the 
ct's transform like conventional Ising gauge fields. This 
is also evident from the structure of the effective action 



obtained by substituting the paramet erization in (4.6) 
and Fig [?] into the effe ctive action in ( 4.1 
terms displayed in (4.1, p~3|) we obtain: 



y); for the 



S3 = / dr 



C2<7iCr 2 0'3a'4A 1 2A i 4 

+C3aia 2 <7 3 a 4: ^i^2 ~ ^1^4 + A*3^4 - ^3^2} 

+C 4 |/il^4 - MlM2 + M2^3 - 

-050-1020-30-4 jl + M1M2M3M4} 

+C6O1CT2O3O4/K1M3 ■ (4.7) 



The terms involving the Ui appear to have the plaque- 
tte form associated with Ising gauge fields. The spatial 
structure of these gauge interactions is made clearer by 
the transformation in Fig |^. Here, we have collapsed 
pairs of sites connected by the diagonal links into single 
sites — we now see that the Oi can viewed as residing on 
the links of a square lattice which is tilted by 45° from 
the original lattice, and their gauge interactions have the 
usual form around elementary plaquettes. 

The fii constitute a separate global Ising degree of free- 
dom associated with the breaking of the symmetry of 
90° spatial rotations between the horizontal and vertical 
directions. In the mean-field theory of the deconfincd 
phase, the state with ^ = 1 corresponds to the state 
with dominant spin correlations at the wavevector (ir, q) 
(say). The degenerate partner state with spin correla- 
tions at (q, 7r) is obtained by the state [ii = (— \) %v , where 
{i x ,i y ) are the Cartesian co-ordinates of the site i. 



So the action S3 describes a Z 2 gauge theory (a) cou- 
pled (rather intricately) to an Ising spin field (//); the 
fi field does not carry a non-zero charge under the Z 2 
gauge transformation. The Z 2 gauge theory can undergo 
a confincment-deconfinement transition (which is related 
by a duality transformation to the magnetic transition 
in an Ising model in three dimensions), corresponding to 
the liberation of spinons upon moving out of the dimer 
phase. In a different sector, the ordering of the fi degrees 
of freedom leads to the appearance of nematic order, and 
the breaking of the symmetry of 90° spatial rotations. 
In the mean-field theory, these two transitions occur at 
the same point i.e. the deconfmcmcnt transition is also 
the point where the spatial rotation symmetry is broken. 
More generally, the interplay between these two poten- 
tially distinct transitions can be addressed by an analysis 
of fluctuations using the action S3 . It does appear pos- 
sible that the two transitions are not simultaneous, and 
that there can be a deconfined phase without any broken 
spatial symmetries; moreover, if there is a simultaneous 
transition in the two sectors, it is likely to be first or- 
der. A more definitive conclusion on these issues must 
await a complete study of the coupled Ising gauge/Ising 
spin theory defined by S3 . We note that these issues con- 
cerning the transition from the confined dimer phase to 
the deconfined helical SRO phase are somewhat different 
from earlier deconfinement transitions^! because here the 
dimer phase does not break any lattice symmetries. 



V. CONCLUSIONS 

The Mott insulator SrCu2 (603)2 is perhaps the only 
example of a spin gap paramagnet on a strongly frus- 
trated two-dimensional lattice (another example of a two- 
dimensional paramagnet is CaA^tOg, but its spin gap is 
realized by dilution and not frustration). To date, it ap- 
pears that the spin gap is realized in a simple decoupled 
dimer groufld state discovered originally by Shastry and 
SutherlandB. Here, we undertook a more detailed study 
of the parameter space of this antiferromagnet, and found 
that other paramagnetic spin gap states are also possi- 
ble. One of these was the plaquette statetj, which appears 
in a region of weaker frustration and commensurate spin 
correlations. The other was a more exotic state with 
"topological order", deconfined 5 = 1/2 excitations, and 
helical spin correlations. The latter state was found to be 
contiguous to the dimer state, and so not too far from the 
physically relevant regime: it appears that SrCu2 (603)2 
is quite close to the boundary of stability of the dimer 
phase. 

Our results suggest exciting possibilities for materials 
obtained by doping SrCu 2 (B03)2 with mobile carriers. It 
is expected that the helical state will be more amenable 
to the motion of charge carriers than the dimer state, 
and so doping may well drive the system into a topologi- 
cally ordered state. Such a state is a prime candidate for 
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superconductivity with the exotic properties associated 
with the proximity of a Mott insulator with deconfincd 
spinons: these include the flux-trapping effect of Santhil 
and Fishera, and a regime of stable hc/e vorticesoE3. 
An experimental effort to dope SrCu2(BC>3)2 (or related 
compounds) therefore appears worthwhile. 
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that the eigenvalues occur in pairs with opposite signs 
(u>^(k), —Lj^(k)) where = 1, • • • ,4, and the matrix M 
has the form 



M 



U -V* 

V u* 



(A4) 



where the U, V are 4x4 matrices associated with the 
positive eigenvalues. The (it, q) LRO phase (xf ^ 0) 
occurs at the wavevector k min — (±7r/2, ±g/2) where the 
eigenenergy vanishes, ie. w(fc m in) = 0. The two linearly 
independent eigenvectors associated with k% = (tt/2, q/2) 
and k 2 — (tt/2, —q/2) can be shown to be 



APPENDIX: 

This appendix will provide a proof of a statement made 
in Section II A 2 on the nature of the saddle point in the 



phase with helical LRO: we will analytically show that in 
the (it, q) and (g, ir) LRO phase the link fields obey the 
following relations: 



Q = P 
\Ri\ = \R 2 \ 



(Al) 



where \Qi\ — Q, |Pj| = P. The reasoning is the same 
in the both (n, q) and (g, n) phase, and we will fix our 
state in the (ir, q) phase for simplicity. In this state, the 
directions of the link fields Qij are shown as in Fig ||. The 
spinon dispersion in this phasencan be obtained from the 
following eigenvalue equationt-3 in momentum space: 



D(k) 



MT 3 ib(k) 

'1 o 
-1 

Al P(k) 



(A2) 



where is the 4x4 unit matrix, uj(k) is a diagonal matrix 
of the bosonic eigenenergies, M is a 8 x 8 matrix whose 
columns are the eigenvectors of the matrix T 3 D(k) and 
the diagonal elements of T 3 t2>(fc) are the corresponding 
eigenvalues, P(k) is a 4 x 4 matrix with the following 
form 



P(k) = 







iJiQsin(k x ) 
-(J 2 i?i/2)e" , ( fe ^ fe ^ 
zJiPsin(fcj,) — 

[J 2 R 1 /2)e l ^- k ^ 
(J 2 P2/2)e J ( fe - +fc ") 

iJiQ sin(fe x ) 

i JiQ sin(fc x ) 



i J\Q sin(fcx) 


iJiP sm(ky) 
{J 2 R 2 /2)e- t ^ +k y^ 

iJ\P s\a(ky) \ 



(A3) 



and A is the Largrange multiplier of the mean-field 
Hamiltonian which we assume to be independent of la-t-i 
tice site i. With this assumption, it can be shownE-3 



(l,ie-^ 2 ,ie-^ 2 ,l,i,- 



-iq/2 



* 2 = (1, «e l « /2 , ie lq/2 , 1, -», e lq/2 , e lq ' 2 , -i)e lk2 ^ (A5) 
respectively. Substituting (or ^2) into (|A^), we have 

A - [,h{P + Q) sin(g/2) + J 2 ^- sin(g)] 



i[J\Q cos(g/2) 



^2^-cos(g)] 
R 2 



[Ji(P + Q) sin(g/2) + J 2 ^ sin(g)] 
R 2 

i[JiQcos(g/2) + J 2 — cos(g)] =0 



(A6) 



We an easily see that R\ = R 2 = R from the above 
conditions. Also, we find that each saddle-point may be 
described by purely real < >. Therefore, we may fix 
the values of A and q in the LRO phase from the above 
condition. 



R 

Ji(P + Q) sin(g/2) + J 2 - sin(g) = A 
R 

JiQcos(g/2) + J 2 —cos(q) = 



(A7) 



To prove P = Q, we need one additional condition from 
the saddle-point equations. The mean-field free energy 
Emf is a function of A, Q, P, R and x a (q) where these 
are independent parameters. The large-N solutions of 
this model are obtained by solving the saddle-point equa- 
tions which set the derivatives of free energy with respect 
to these independent variables to be zero. Notice that q is 
also an independent parameter. The additional condition 
we need comes from the saddle-point equation associated 
with q. It is given by 



dE 



MF 



0. 



(A8) 



The only explicit q dependence in the free energy is in 
the bose condensate variables x a (q). This piece of free 
energy is given byl3 

E M = E If [-Qij^a>xt(q)xf (q) + H.c] , (A9) 

i>j 
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where e CTCT ' is the SU(2) antisymmetric e tensor, and 
a, a 1 =|, |. The condensates xf(q) must be the linear 
combinations of the eigenvectors and \&2 associated 
with the zero mode: this introduces two complex num- 
bers ci, C2, with only the value of \c\ 2 + 1 C2 1 2 fixed by the 
saddle-point equationsLJ. Working out the orientation of 
the condensate at every lattice site over the unit cell, the 
condensates can be written in the form 



ci 



B - 



C - 



C2 

- Cl e- l i/ 2 
-ic* 2 e~ lq / 2 

-icU~ iq - 



c 2 e lq / 2 
ic1e lq/2 



1 

c 2 e lq 



c 2 e 



r,* p -iq/2 



iq/2 
,iq/2 



(A10) 



By substituting ( A10| ) into (A.9), we can explicitly work 
out E x (q) . It is given by 



E x (q) = -[Ji(P + Q)sm{q/2) + 

Now the saddle-point condition 

dE MF t P + Q , , 
— = J 1 ^cos(,/2) + 



R 

J2^sm(q)}(\ Cl \ 2 



jA§| ) becomes 
R 

J 2 —cos(q) = 0. 



+ M 2 ). 

(All) 
(A12) 



Combining (A7) and (A12), we have P = Q 
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FIG. 2. The four sites of the unit cell (labeled A, B, C and 
D), and the 10 link variables Qij. 
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J2 / (Jl + J2) 

FIG. 3. Large N phase diagram of the Sp(N) Shas- 
try-Sutherland model, Eq. [l.l| , as a function of J%/(J\ + J2) 
and 1/5. The five phases are described in Section [n]. The 
LRO phases break spin-rotation symmetry: the spin order 
is collinear and commensurate in the (71", n) LRO phase, and 
helical and incommensurate in the (n, q) LRO phase. The 
SRO phases preserve spin rotation invariance. In the (n, n) 
SRO only the horizontal and vertical Qij are non-zero in the 
large N theory-fluctuations lead to broken translational sym- 
metry in one of the states shown in Fig |(| (A state with 
co-existing47r, tt) LRO and plaquette order is also allowed by 
the theoryt-j beyond the large N limit (not shown above), and 
there is evidence-that this occurs on a frustrated square lattice 
antiferromagnetc3.) The dimer phase has only the diagonal 
Qij non-zero in the large TV theory. The (n,q) SRO phase 
has all the Qij non-zero: this phase has topological order and 
deconfined spinons. 
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FIG. 4. Three states of the Hilbert space of the quantum 
dirtier model. There are off-diagonal matrix elements in the 
effective Hamiltonian which connect state (a) to state (b), 
and state (a) to state (c), by a resonance between pairs of 
horizontal and vertical dimers around a plaquette. The latter 
matrix element differs from the former because only the lat- 
ter has a diagonal link across the resonating plaquette. Also 
shown are the corresponding values of the heights, h a , on the 
sites of the dual lattice. 
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FIG. 6. Phase diagram of (3.6) as a function of the pa- 
rameters A and B; this model describes fluctuations in the 
(-7T, 7r) SRO phase of Fig ^. The thick line is a first order tran- 
sition, while the thin lines are second order. The plaquette 
and spin-Peierls states are shown, with the different line-styles 
representing distinct values of (Si • Sj) across the links. 
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FIG. 5. The four dual sublattices upon which the height 
offsets take the values Cyv — 0, £x = 1/4, £y = 1/2, and 
Cz = 3/4. 
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FIG. 7. A section of the Shastry-Sutherland lattice. We 
have labeled sites around the central plaquette to enable the 
discussion in Section ^ of the various terms in the Zi gauge 
theory of the transition from the dimer state to the (n, q) SRO 
phase with spinon deconfinement. 
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FIG. 8. A deformation of the Shastry-Sutherland lattice 
which exposes the structure of the Z2 gauge theory. Pairs 
of sites across a diagonal bond have been compressed into 
a single site. Four of the sites carry pairs of sites labels, 
corresponding to the original site numbers in Fig |^. The Z2 
Ising gauge fields on some of the links are indicated, with a 
notation corresponding to the degrees of freedom in Fig Fl 
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